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Recently, NASA, FAA, and other organizations have focused their 
attention upon the possible effects of rain on airfoil performance. Rhode 1 carried 
out early experiments and concluded that the rain impacting the aircraft 
increased the drag. Bergrum2 made numerical calculation for the rain effects on 
airfoil. He claimed that thin airfoils of different sections having the same exposed 
frontal areas, will have approximately the same rates of water-drop impingement 
at high speeds, but the distribution of water-drop impingement will be different 
for each section. Luers and Haines 3 did analytic investigation and found that 
heavy rain induces severe aerodynamic penalties including both momentum 
penalty due to the impact of the rain and a drag and lift penalty due to rain 
roughening of the airfoil and fuselage. More recently, Hansman and Barsotti* 
performed experiments and declared that performance degradation of an airfoil 
in heavy rain is due to the effective roughening of the surface by the water layer. 
Hansman and Craig 5 did further experimental research at low Reynolds number. 
They concluded that the initial effect of rain is to cause premature boundary layer 
transition near the leading edge. 

E. Dunham 5 made critical review for the potential influence of rain on 
airfoil performance. Dunham 7 et al. carried out experiments for the transport type 
airfoil and concluded that there is a reduction of maximum lift capability with 
increase in drag. There is a scarcity of published literature in analytic research 
of two-phase boundary layer around an airfoil. Although Henry 8 et al. presented 
a technical paper entitled "A Von Karman Integral Approach to a Two-Phase 
Boundary Layer", yet their main assumption that the existence of zero shear 
stress at the wall and on the interface is not physically realistic. Most recently 
Bilanin 9 attempted an analytic investigation. 

He assumed that the ejacta layer thickness is constant in his preliminary 
report. This assumption is quite doubtful. The present author attempts to 
improve the analytic research. The following assumptions are made: 

1. The fluid flow is non-steady, viscous, and incompressible. 

2. The airfoil is represented by a two-dimensional flat plate. 

3. There is only laminary boundary layer throughout the flow region. 

Under the usual boundary layer approximation, there obtains, 
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Under the usual boundary layer approximation, there obtains, 


For the liquid: dii_ + dx = a. Yi sinf) (1) 

Eq. of continuity dx dy 8 

Eq. of momentum Ou + «au + \> du = 1 d 2 u +a Vj^sinpcosP (2) 

at ax ay (Rj t a y z s 

Where (R e ) e = PwVrL = Reynolds number of the liquid water 

M- 

a = Wl = ration of liquid water content to water density 

P w 
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For the fog: 


Equation of continuity 

shL + 
dx 

11 

o 



(3) 

Equation of momentum 
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The initial and boundary conditions are: 


at t = 0 'll = U(x,0) V=0 

for the liquid phase, at y=0, u=v=0 (No slip condition) 

at the interface y=8 

= (% 

. (\k = (V)/ x 

as y — > 'll = U (x,x) 

The above set of partial differential equations is non-linear in nature. An exact solution 
is not possible. The said set of partial differential equations is transformed into a set of 
finite difference equations. Using fortran language. An numerical solution is expected. 
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